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The v e flux and spectrum have been measured at a distance of about 800 m from the reactors of 
the Palo Verde Nuclear Generating Station using a segmented Gd-loaded liquid scintillator detector. 
Correlated positron-neutron events from the reaction v e p— >e + n were recorded for a period of 200 d 
including 55 d with one of the three reactors off for refueling. Backgrounds were accounted for by 
making use of the reactor-on and reactor-off cycles, and also with a novel technique based on the 
difference between signal and background under reversal of the e + and n portions of the events. 
A detailed description of the detector calibration, background subtraction, and data analysis is 
presented here. Results from the experiment show no evidence for neutrino oscillations. v c — > P x 
oscillations were excluded at 90% CL for Am 2 > 1.12 x 10" 3 eV 2 for full mixing, and sin 29 > 
0.21 for large Am 2 . These results support the conclusion that the observed atmospheric neutrino 
oscillations does not involve v e . 

PACS 13.15.+g, 14.60.Lm, 14.60.Pq 



I. INTRODUCTION 

Results of a long baseline study of D c oscillations at 
the Palo Verde Nuclear Generating Station are reported 
here. The work was motivated by the observation of an 
anomalous atmospheric neutrino ratio v^jv^ reported in 
several independent experiments that can be inter- 
preted as v^-v c oscillations requiring large mixing. The 
mass parameter suggested by this anomaly is in the range 
of 1CP 2 < Am 2 < 10~ 3 eV 2 for two flavor neutrino oscil- 
lations. 

The quantity Am 2 , defined as the difference between 
the square of the masses of the mass eigenstates, and the 
mixing parameter 6 are related to the transition probabil- 
ity P for two-flavor v a — > i/f, oscillations (see, for example, 
1) by: 

„ / \ .2^.3 /L27Ato 2 L\ 

P osc {v a -» vb) = sin 2 20sin 2 ( J , (1) 

where E v (MeV) is the neutrino energy, L (m) is the 
source-detector distance, and Am 2 is measured in eV 2 . 

Exploring Am 2 down to 1CP 3 eV 2 requires that the 
quantity L/E v (m/MeV) has a value of around 200. For 
reactor neutrinos (E v ^5 MeV), a baseline of L ~1 km is 
adequate. Reactor experiments are generally well suited 
to study v e oscillations at small Am 2 ; however, they are 
restricted to the disappearance channel v e — > v x . 

Reactor antineutrinos have been used for oscillation 
studies with ever increasing Am 2 sensitivity since 1981 
UK . All of the experiments are based on the large cross 
section inverse beta decay reaction, ^ e p— >e + n. The corre- 
lated signature, a positron followed by a neutron capture, 



allows significant suppression of backgrounds. As the re- 
actor D e yield and spectra are well known a "near 
detector" is not required. It is, however, important to 
control well the detector efficiency and backgrounds. 

The mentioned considerations have led to the design 
of the Palo Verde and Chooz JjJ experiments, which have 
similar Am 2 sensitivities. While both experiments have 
pursued their goal of exploring the unknown region of 
small Am 2 , recent data from Super-Kamiokande || fa- 
vor the Vp — > v x oscillation channel over — > v c . This 
paper reports in greater detail results presented earlier || 
and describes the detector calibration, background sub- 
traction, and data analysis techniques used to extract 
results on neutrino oscillations. 

II. THE EXPERIMENT 

A. The detector 

The Palo Verde Nuclear Generating Station in Ari- 
zona, the largest nuclear power plant in the USA, con- 
sists of three identical pressurized water reactors with 
a total thermal power of 11.63 GW. The detector is lo- 
cated at a distance of 890 m from two of the reactors 
and 750 m from the third at a shallow underground site. 
The 32 meter-water-equivalent overburden entirely elim- 
inates any hadronic component of cosmic radiation while 
reducing the cosmic muon flux to 22 m s . In or- 
der to reduce the ambient 7-ray flux in the laboratory 
all materials in and surrounding the detector were se- 
lected for low activity. The laboratory walls were built 
with an aggregate of crushed marble, selected for its low 
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content of natural radioisotopes. Concentrations of 170, 
750, and 560 ppb for 40 K, 232 Th, and 238 U were measured 
in the concrete resulting in a tenfold reduction of 7-ray 
flux when compared with locally available aggregate. A 
low 222 Rn concentration of about 20 Bq/m 3 in the lab 
air was maintained with forced ventilation. Temperature 
and humidity were controlled to ensure stable detector 
operation. 

The segmented detector, shown in Fig. |l|, con- 
sists of a 6x11 array of acrylic cells dimensioned at 
900 cmxl2.7 cmx25.4 cm and filled with a total of 
11.34 tons of liquid scintillator. A 0.8 m long oil buffer 
at the ends of each cell shields the central detector from 
radioactivity originating in the photomultiplier tubes 
(PMTs) and laboratory walls. The cells were made by 
cutting and bonding large 0.62 cm thick acrylic sheets. 
The total acrylic mass in the detector is 3.48 tons. Each 
cell is individually wrapped in 0.13 mm thick Cu foil to 
ensure light-tightness and is viewed by two 5-inch low 
activity PMTs |[o]], one at each end, housed in mu-metal 
boxes. The target cells are suspended on rollers held in 
place by thin sheet metal hangers. All structural mate- 
rials were dimensioned as lightly as possible to minimize 
dead material between cells. Each cell can be individu- 
ally removed from the mechanical structure for mainte- 
nance. The detector is oriented such that the v e flux is 
perpendicular to the long axis of the cells. 

The liquid scintillator is composed of 36% pseudoc- 
umene, 60% mineral oil, and 4% alcohol, and is loaded 
with 0.1% Gd by weight. This formulation was chosen 
to yield long light transmission length (11.5 ± 0.1 m at 
440 nm), good stability, high light output, and long term 
compatibility with acrylic. Details of the scintillator de- 
velopment have been published elsewhere pd| ] . 

The central volume is surrounded on the sides by a 1 m 
buffer of high purity deionized water (about 105 tons) 
contained in steel tanks which, together with the oil 
buffers at the ends of the cells, serve to attenuate gamma 
radiation from the laboratory walls as well as neutrons 
produced by cosmic muons passing outside of the de- 
tector. The low Z of water minimizes the neutron pro- 
duction by nuclear capture of stopped muons inside the 
detector and has a high efficiency for neutron thermal- 
ization. 

The outermost layer of the detector is an active muon 
veto counter, providing An coverage. It consists of 32 
twelve meter-long PVC tanks (from the MACRO exper- 
iment |l2|]) surrounding the detector longitudinally, and 
two endcaps. The endcaps are mounted on a rail system 
to allow access to the central detector. The horizontal 
tanks are read out by two 5-inch PMTs at each end; the 
vertical tanks are equipped with one 8-inch PMT at each 
end while the endcaps use 3-inch PMTs. The liquid scin- 
tillator used in the veto is a mixture of 2% pseudocumene 
and 98% mineral oil, with a light attenuation length at 
440 nm in excess of 12 m. 
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FIG. 1. The Palo Verde detector. One of the 66 target 
cells with PMTs, calibration LEDs, and optical fiber flashers 
is shown lengthwise at the bottom. 



A schematic of the central detector's front-end elec- 
tronics is shown in Fig. |[ Each channel can be digi- 
tized by either of two identical banks of electronics. The 
dual bank system allows both parts of the sequential in- 
verse beta decay event to be recorded with no deadtime 
by switching between banks. Due to the large dynamic 
range of energy in the data of interest (40 keV to 10 
MeV, or 1 to 250 photoelectrons typically), each PMT 
has both a dynode and anode output connected to ADCs, 
as well as three discriminator thresholds for the trigger 
and TDCs. The higher TDC threshold serves to avoid 
crosstalk from large signals in adjacent channels while 
the lower threshold allows timing information to still be 
available at the single photoelectron level. The relative 
time of arrival from each end of a cell is used to recon- 
struct longitudinal position. The measured PMT pulse 
charge at each end, corrected for light attenuation based 
on the distance traveled in the cell, allows energy recon- 
struction. 

Each cell is connected to the trigger via the OR of the 
discriminated signals from the two PMTs. Signals are 
tagged according to two thresholds: a high threshold cor- 
responding to ^600 keV for energy deposits in the middle 
of the cell and a low threshold corresponding to ~40 keV, 
or one photoelectron at the PMT. The low trigger thresh- 
old also serves as the lower TDC threshold. The trigger, 
which has a decision time for each event of around 40 ns, 
uses a Field Programmable Gate Array to search for pat- 
terns of energy deposits in the central detector, and can 
be reprogrammed easily to change trigger conditions as 
needed for calibrations Jl3|. 
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FIG. 2. Schematic of the front-end electronics. The charge 
and timing of PMT pulses are read out by two banks of ADCs 
and TDCs. The timing information is discriminated with 
two thresholds to avoid crosstalk and retain dynamic range. 
Trigger discriminator signals from the two PMTs in each cell 
are OR'ed and input to the trigger circuit. 



A veto signal disables the central detector trigger for 
10 (is following the passage of a muon to avoid most 
related activity. Typical veto rates are ~2 kHz. With 
each event, the time and hit pattern of the previous muon 
in the veto counter is recorded along with information as 
to whether or not the muon passed through the target 
cells. The veto inefficiency was measured to be (4±1)% 
for stopping muons (one hit missed) and (0.07±0.02)% 
for through-going muons (two hits missed). We note that 
the small size of this second quantity with respect to the 
first is due to correlations between incoming and outgoing 
muons as confirmed by a simple Monte Carlo model. 



B. The v c signal 

The D c signal is detected via the reaction v e p— >ne + as 
illustrated in Fig. [I?] further below along with the domi- 
nant backgrounds. Signal events consist of a pair of time- 
correlated subevents: (1) the positron kinetic energy ion- 
ization and two annihilation 7's forming the prompt part 
and (2) the subsequent capture of the thermalized neu- 
tron on Gd forming the delayed part. By loading the 
scintillator with 0.1% Gd, which has a high thermal neu- 
tron capture cross section, the neutron capture time is 
reduced to ~27 /^s from ~170 /is for the unloaded scintil- 
lator. Furthermore, Gd de-excites by releasing an 8 MeV 
7 cascade, whose summed energy gives a robust event 
tag well above natural radioactivity. In contrast, neu- 
tron capture on protons releases only a single 2.2 MeV 
7- 

Background is rejected at trigger level using the detec- 
tor segmentation by looking for coincidences of energy de- 
posits matching the pattern of inverse beta decay. Each 



of the subevents of a v e signal is triggered by scanning the 
detector for a pattern of three simultaneous hits in any 
3x5 subset of the cell array. This threefold coincidence, 
called a triple, must consist of at least one high trigger hit, 
due to either the positron ionization or neutron capture 
cascade core, and at least two additional low trigger hits, 
resulting from either positron annihilation 7's or neutron 
capture shower tails. The use of identical trigger require- 
ments for the two triples is found to give rise to close to 
an optimal signal to noise ratio. Five fj,s after finding an 
initial triple, the trigger begins searching for a delayed 
triple. The blank time suppresses possible false signals 
from PMT afterpulsing. If two triples are found within 
450 fj,s of each other, the candidate v c event is digitized 
for offline analysis. 



C. Expected v c interaction rate 

In order to calculate the expected v c interaction rate in 
the detector, the status of the three reactors is tracked 
daily, and the fission rates in the cores are calculated 
based on a simulation code provided by the manufacturer 
of the reactors. This code uses as input the power level of 
the reactors, various parameters measured in the primary 
cooling loop, and the original composition of the core fuel 
elements. 

The output of the core simulation has been checked 
by measuring isotopic abundances in expended fuel el- 
ements in the core; errors in fuel exposure and isotopic 
abundances are estimated to cause < 0.3% uncertainty in 
the i> e flux estimate. Of the four isotopes — 239 Pu, 241 Pu, 
235 U, and 238 TJ — whose fissions produce virtually all of 
the thermal power as well as neutrinos, measurements of 
the neutrino yield per fission and energy spectra exist for 
the first three [ p^Jl^] . The 238 U yield, which contributes 
11% to the final v e rate, is calculated from theory p6[ . 
When the same theoretical method was used to calculate 
the spectrum from the other three isotopes, the theory 
agreed with experimental results within 10%. The con- 
tribution of 238 U fission to the overall uncertainty in D c 
rate is therefore expected to be ~1%. 

This calculated v e flux is then used to compute the 
expected rate of v e candidates at the detector as a 
function of the oscillation parameters Am 2 and sin 2 29: 



n p / dEc c a(E PB )r](Ep B )x 



* T Pe ,j(g g J (1 - P OSC;i (Am 2 , sin 2 29, U, E Bc )) 

i=l 



4tt£ 2 



(2) 



where a^E^) is the inverse beta decay cross section [ p^f] , 
rj(Ev m ) is the (energy dependent) detector efficiency, n p is 
the number of target free protons, and 2p c .i is the source 
strength of reactor i at distance Li with oscillation proba- 
bility P sc,i- In Fig. U we show the energy spectrum of the 
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FIG. 3. v,, flux, inverse beta decay cross section, and i> e 
interaction spectrum at the detector. 



u e 's emitted by a reactor, the v c (energy) differential cross 
section in the detector and the actual interaction rate in 
the detector target before detector efficiency corrections, 
referred to here as R Pc (obtained by setting ry(_E p J = l). 
The energy spectrum actually measured in the detector 
is the energy of the positron created by the inverse beta 
decay. This spectrum is approximately E„ e — 1.8 MeV, 
slightly modified by the kinetic energy carried away by 
the neutron (~50 keV). 

Previous short baseline experiments which measured 
the rate of i> e emission by reactors have found good 
agreement between calculated and observed neutrino 
flux by using largely the same method of calcula- 
tion. A high statistics measurement at Bugey ||, in 
particular, found excellent agreement both in spectral 
shape (x 2 /n.d.f.=9.23/ll) and in absolute neutrino yield 
(agreement better than 3%, dominated by systematic er- 
rors). These previous generation experiments prove that 
the reactor antineutrino spectrum, i.e. the v B flux at the 
distance L—Q, is well understood. 

The expected v c interaction rate in the whole target, 
both scintillator and the acrylic cells, is plotted in Fig. |] 
for the case of no oscillations from July 1998 to October 
1999. Around 220 interactions per day are expected with 
all three units at full power. The periods of sharply re- 
duced rate occurred when one of the three reactors was 
off for refueling, the more distant reactors each contribut- 
ing approximately 30% of the rate and the closer reactor 
the remaining 40%. The short spikes of decreased rate 
are due to short reactor outages, usually less than a day. 
The gradual decline in rate between refuelings is caused 
by fuel burnup, which changes the fuel composition in the 



core and the relative fission rates of the isotopes, thereby 
affecting slightly the spectral shape of the emitted v c flux. 
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FIG. 4. The calculated v a interaction rate in the detector 
target. The two long periods of reduced flux from reactor refu- 
elings were used for background subtraction. The decreasing 
rate during the full power operation is a result of the changing 
core composition as the reactor fuel is burned. 



III. CALIBRATION 

In order to maintain constant data quality during run- 
ning, a program of continuous calibration and monitor- 
ing of all central detector cells is followed. Blue LEDs 
installed inside each cell are used for relative timing and 
position calibration. Optical fibers at the end of each 
cell, also illuminated by blue LEDs, provide information 
about PMT linearity and short term gain changes. LED 
and fiberoptic scans are performed once a week. Ra- 
dioactive sources are used to map the light attenuation 
in each cell, for absolute energy calibration, and to de- 
termine detection efficiencies for positrons and neutrons. 
A complete source scan is undertaken every 2-3 months. 



A. LED and optical fiber calibrations 

As seen in Fig. [I], every cell of the central detector 
has two LEDs, one at each end at a distance of 90 cm 
from the PMTs. These blue LEDs, which provide fast 
light pulses with a rise time comparable to scintillation 
light, are used for timing calibrations needed for position 
reconstruction along the cell's axis. 
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The difference in pulse arrival time between the two 
PMTs of a cell At is described as a function of the po- 
sition z with an effective speed of light c e s, an offset zq 
and a small nonlinear correction /(1/Q): 



At = (z-z )/c efl + f{l/Q n ,l/Q f ). 



(3) 



The correction /(1/Q n , 1/Q/), a function of both near 
and far PMT pulse charge Q, describes the dependence 
of the pulse height due to time-walk in the leading edge 
discriminators used in the front-end electronics. To ex- 
tract these calibration parameters and compensate for 
the time-walk effect, a third order polynomial is ht to 
At versus 1/Q (see Fig. ||). The intercepts at 1/Q = 
for the two LED positions provide c e g and z , while the 
slopes are used to parameterize the time-walk correction. 
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FIG. 5. Time difference At as function of 1/Q for near- and 
far-end LEDs of a cell, illustrating the time-walk effect. For 
values of Q close to the discriminator threshold the time-walk 
correction is kept at a constant (maximum) level. 

In order to check the suitability of longer wave- 
length 470 nm LED light to measure timing proper- 
ties of ~425 nm scintillation light, data taken with a 
228 Th source at several longitudinal positions were re- 
constructed with the LED timing calibration parameters. 
Comparing the reconstructed positions with the actual 
source positions, the effective speed of light measured 
with the LED system was found to be on average 3.6% 
lower than that with the sources. A simulation of the 
light transport in a cell with various indices of refraction 
and attenuation lengths of the scintillator suggested that 
the small discrepancy in At between LED and scintilla- 
tion light was due to the difference in attenuation length. 
The correction factor was found to be constant over sev- 
eral months. Weekly LED scans are therefore used to 



correct for short term variations in At and a constant 
correction factor is applied to the effective speed of light. 

The fiberoptic system includes 15 blue LEDs, each illu- 
minating a bundle of 12 fibers. The light output of each 
LED is measured in two independent reference cells with 
PMTs checked to be linear over the whole dynamic range 
of the LEDs. By taking a run which scans through all 
light intensities and mapping each PMT's response rela- 
tive to the reference cells, the nonlinear energy response 
of the PMTs is calibrated. Low intensities are used to 
determine the single photoelectron gain of each PMT, 
which is used to correct for changes from the nominal 
gain setting of 4 x 10 7 . 



B. Scintillator transparency and energy scale 
calibration 



In addition to weekly LED and fiberoptic calibrations, 
the energy response of the scintillator is measured every 
three months using a set of sealed radioactive sources. 
Eighteen 2.4 mm diameter tubes run along the length 
of the detector allowing insertion of the sources adjacent 
to any cell at any longitudinal position. The response 
of each PMT as a function of longitudinal position is 
measured by recording the Compton spectrum from the 
2.614 MeV 7 of a 228 Th source at seven different locations 
along each cell. 

Monte Carlo simulation found that the half maximum 
of a Gaussian function fitted to the Compton spectrum is 
relatively independent of resolution; this point is there- 
fore used as the benchmark of the cell response. The 
response versus distance from the PMT, shown in Fig. ^ 
for one cell, is then fit to the phenomenological function 
e-xjp(pQ+piz)+exp(p2+P3z)/ z, where z is source longitu- 
dinal distance from the PMT. The effective attenuation 
length of the scintillator (including multiple total reflec- 
tion on the acrylic walls) is generally between 3-4 m and 
over a year was found to change on average ~1 mm/day, 
demonstrating that the Gd scintillator was remarkably 
stable. 

The overall energy scale was determined from the po- 
sition of the 1.275 MeV peak of a 22 Na source, and then 
verified by taking data with several 7 sources in differ- 
ent energy ranges: 137 Cs (0.662 MeV), 65 Zn (1.351 MeV), 
228 Th (2.614 MeV), and the capture of neutrons (8 MeV) 
from an Am-Be source. The gamma cascade from neu- 
tron capture was modeled according to measurements 
of the emitted spectrum |TJ|]. In contrast to homoge- 
neous detectors which measure total absorption energy 
peaks, 25% of the detector target mass consists of the 
inert acrylic of the cell walls, which absorbs some en- 
ergy. The Monte Carlo simulation was therefore used 
to find the correct final distributions of energy detected 
from single and multiple scattering of the 7's. The to- 
tal energy reconstructed for data and Monte Carlo for 
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FIG. 6. Effect of aging on Gd loaded scintillator. Top: 
228 Th Compton edge position at seven different longitudinal 
locations along a typical cell. The two curves are from cali- 
brations taken a year apart. The curves are normalized at the 
location nearest to the PMT. Bottom: Effective attenuation 
lengths for all 66 cells from the two calibrations. 



each source is plotted in Fig. ^ The data were matched 
with Monte Carlo simulation for the 22 Na spectrum in 
Fig. |l2] to find the overall energy scale and to the spectra 
in Fig. ^ to assure that the scintillator response is linear 
over the energies of interest. The light yield after PMT 
quantum efficiency was found to be ~50 photoelectrons 
per MeV in the center of the cells. The agreement for 
three of the four sources in Fig. fj] is good, the excep- 
tion being 228 Th, in which the data has a consistently 
higher Compton scattering peak than Monte Carlo pre- 
dicts. This discrepancy is consistent across all the data 
taken and therefore does not affect the scintillator trans- 
parency calibration. 



C. Monte Carlo simulation 

The v c efficiency of the detector is a relatively strong 
function of event location in the detector and, to a lesser 
extent, of time due to scintillator aging. A further com- 
plication comes from the trigger efficiency being a func- 
tion of threshold (voltage) while only energy (charge) is 
measured. For this reason a Monte Carlo model which 
included a detailed simulation of the detector response, 
including the PMT pulse shape, is used for an estimate of 
the overall efficiency for D c detection. A variety of mea- 
surements was performed to crosscheck that the Monte 
Carlo accurately models the detector response. 
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FIG. 7. The total energy reconstructed for various sources 
compared for data and Monte Carlo. 



The physics simulation program is based on GEANT 
3.21 p9|| . This code contains the whole detector geome- 
try and simulated the energy, time, and position of en- 
ergy deposits in the detector. Hadronic interactions are 
simulated by GFLUKA ]20| and the low energy neutron 
transport by GCALOR Scintillator light quenching, 

parameterized as a function of ionization density, is in- 
cluded in the simulation p2[ . 

The event reconstruction program reads the output of 
this physics simulation and then applies the second step 
of the Monte Carlo, the simulation of the detector re- 
sponse as PMT pulses which are then converted into time 
and amplitude digitizations and trigger hits. A logical 
scheme of this detailed detector simulation is shown in 
Fig. |. 

The calibrations discussed above empirically provide 
the scintillator light yield (photoelectron/MeV) and at- 
tenuation function for each cell, which in turn provide 
the number of photoelectrons iV p e . expected for a simu- 
lated energy deposit. The total charge of the pulse then 
follows from sampling a Poisson distribution with mean 
A^p.e. and folding the number of simulated photoelectrons 
Np. e . with the PMT's nominal gain, first stage gain vari- 
ance (10iVp. e .), and cell-to-cell energy scale calibration 
uncertainty (10%). 

To simulate the pulse shape, an arrival time is as- 
signed to each photoclcctron, and individual photoclcc- 
tron pulses (whose shape is derived from real data) are 
summed into a final pulse. The calculated arrival time of 
each photon is a combination of two processes, scintillator 
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FIG. 8. Schematic of the Monte Carlo detector simulation 
program. In order to convert energy deposited into accurate 
threshold simulation, calibration information from the detec- 
tor was used to reconstruct PMT pulses. 



de-excitation and propagation along the cell. The latter 
distribution is parameterized by the distance traveled to 
the PMT, larger distances giving larger variances, using 
a light transport simulation of 2 x 10 7 photons. The re- 
sulting pulse is then analyzed to extract TDC and trigger 
hits. 

The Monte Carlo threshold simulation, position recon- 
struction, and positron and neutron efficiency predictions 
were checked using calibration data. The trigger thresh- 
old simulation for each cell was compared to data taken 
with a 22 Na (3 + source near the center of each cell. The 
trigger conditions were loosened for these data, a single 
low hit producing a trigger and the event tagged if a high 
threshold was crossed. By plotting the reconstructed en- 
ergy for each event versus the efficiency for a high trigger 
tag, an effective high trigger threshold in MeV for that lo- 
cation in the cell was determined. The low threshold was 
measured similarly. The Monte Carlo pulse shape param- 
eters were tuned to these data. A typical cell's trigger 
threshold efficiency as a function of energy is shown in 
Fig. ^ for both data and the Monte Carlo. The trigger 
threshold, defined as the energy at 50% efficiency, is also 
plotted for all 66 cells. On average, the thresholds agree 
to within 1%. 

TDC thresholds were checked by the same algorithm, 
plotting the threshold hit efficiency versus reconstructed 
energy. A more direct check of the TDC simulation, how- 
ever, compares the position reconstruction for data and 
Monte Carlo simulation. Fig. [To] shows the longitudi- 
nal position of the third largest energy deposit in each 
event for a 22 Na calibration run, representing the po- 
sition reconstruction of the energy deposited by one of 
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FIG. 9. A comparison of the trigger thresholds from data 
and Monte Carlo. The data were taken with a 22 Na source at 
the center of each cell. The top portion shows the efficiency of 
the trigger threshold for a typical cell as a function of energy 
deposited; the bottom shows the energy of the 50% efficiency 
threshold for all 66 cells. 



the two positron annihilation 7's. Since these energy de- 
posits tend to be small (~100 keV), some fraction of them 
have one or both PMT's responses below the TDC low 
threshold. These events constitute the tails of the distri- 
bution in Fig. |l0] since only the relative signal amplitude 
was used for position reconstruction. The narrower cen- 
tral peak is populated by events with TDC information 
available. The simulation and data agree well, in both 
resolution and relative frequency of the two cases. 



D. v c detection efficiency 

The absolute efficiency of the detector for positron an- 
nihilations and neutron captures was verified using 22 Na 
and Am-Be sources respectively. The 22 Na source emits 
a 1.275 MeV primary 7 which is accompanied 90% of the 
time by a low energy positron which annihilates in the 
source capsule. The primary 7 can mimic the positron 
ionization of a low energy 9 e event. This deposit, in con- 
junction with the positron's annihilation 7's, closely ap- 
proximates the positron portion of a D c event near the 
trigger threshold. 

In two rounds of data taking, 10 months apart, the 
22 Na source was inserted into the central detector at 35 
locations chosen to provide a sampling of various dis- 
tances from the PMTs and edges of the fiducial volume. 
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FIG. 10. A comparison of the reconstructed Z position of 
an annihilation gamma energy deposit from a 22 Na run, be- 
tween data and Monte Carlo. The long tails are due to lower 
energy deposits which do not have TDC for position recon- 
struction, in which case relative signal amplitudes are used, 
resulting in worse resolution. Statistical errors in the plot are 
negligible. 
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FIG. 11. Comparison of data and Monte Carlo detection 
efficiency for Am-Be and 22 Na source runs at various loca- 
tions. Locations at the edges of the detector tend to have 
lower efficiencies. 



The source activity is known to 1.5%, allowing determi- 
nation of an absolute efficiency. After applying the of- 
fline selections used for v c prompt triples and correcting 
for detector DAQ deadtime, the measured absolute effi- 
ciency was compared with the Monte Carlo prediction; 
the results are summarized in the top portion of Fig. [ll]. 
Good agreement is seen in the average efficiency over all 
runs, and run by run agreement was 11%. 

The energy spectra predicted by the simulation and 
measured in the data for the 22 Na runs were compared. 
The total energy seen in all cells and the energy detected 
in the three most energetic hits are plotted in Fig. |lj. 
The trigger thresholds can be seen in the spectra: the 
high trigger threshold is the rising edge at around 0.5 
MeV in the spectrum of the most energetic hit (Ei), and 
the low trigger threshold is the rising edge at around 50 
keV of the third most energetic hit (£3). 

A similar procedure was used to check the neutron cap- 
ture detection efficiency. The Am-Be neutron source is 
attached to one end of a thin (7.5mm) NaI(Tl)-detector, 
which tagged the 4.4 MeV 7 emitted in coincidence with 
a neutron. The Nal(Tl) tag forces the digitization of the 
4.4 MeV 7 as the prompt part of an event and opens 
a 450 (is window for neutron capture; this is the same 
coincidence window used in the v c runs. 

All neutron cuts used for the 9 e data selection were ap- 
plied, and the resulting detection efficiency was corrected 




E 2 (MeV) E 3 (MeV) 

FIG. 12. The Monte Carlo simulation and the data com- 
pared for the spectra of total energy and first, second, and 
third most energetic hit (Stotai, E\, E2, and E3) digitized in 
the 22 Na calibration runs. 

for detector deadtime and a small random coincidence 
background. On average, the Monte Carlo efficiency pre- 
dictions agrees well over the 25 locations tested with a 
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run by run agreement of better than 4%, as shown in the 
bottom of Fig. O. 

As with the 22 Na runs, the energy spectra predicted by 
the simulation and measured in the data were compared. 
The total energy seen in all cells and the energy detected 
in the three most energetic hits is plotted in Fig. [l^. 
Note the small peak in £?total at ^2 MeV arising from 
neutrons being captured on hydrogen. The differences 
in data versus Monte Carlo spectra for 22 Na and Am-Bc 
were taken into account in estimating systematic errors. 
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FIG. 13. The Monte Carlo simulation and the data com- 
pared for the spectra of total energy and first, second, and 
third most energetic hit (-Etotai, E\, E2, and £3) digitized in 
the Am-Be calibration runs. 

The Am-Be source emits neutrons with kinetic ener- 
gies up to 10 MeV, creating proton recoils in the detector 
scintillator in coincidence with the Nal(Tl) induced trig- 
ger. By digitizing any energy deposits seen during the 
neutron release, the high ionization density of these re- 
coiling protons was used in setting the parameters which 
control scintillator light quenching in the simulation. 

The above crosschecks verify our ability to accurately 
generate the events, model the detector response, recon- 
struct the events, and correctly calculate the livetime 
of the data acquisition (DAQ) system. Taken together 
these procedures complete the task of estimating our v e 
efficiency. 

The Monte Carlo simulation for D c events models the 
expected interactions throughout the entire target, in- 
cluding the acrylic walls of the cells, since there is signif- 
icant efficiency for inverse beta decay originating in the 
acrylic. The Monte Carlo simulation yields an average 
efficiency over the entire detector as a function of u e en- 



ergy. The efficiency from the simulation is folded with 
the incident z/ e spectrum (which may be distorted by os- 
cillations depending on the hypothesis tested), to get the 
effective efficiency. 

E. An independent reconstruction and Monte Carlo 

A parallel and independent event reconstruction and 
the simulation of the detector response has been de- 
veloped. This second version follows the same general 
outline of detector calibration, event reconstruction, and 
simulation described above, but differs in the algorithms 
and parameterizations used. Major differences include: 

• The functional form for the scintillator light attenu- 
ation is the sum of two exponentials exp(po +p\z) + 
exp(p2 +P3z), without z^ 1 in the second term. 

• The cell response benchmark is the 70% maximum 
rather than half maximum of the fitted Compton 
scattering spectrum. 

• A different parameterization is used for the linear- 
ity correction of the dynode signals. 

• The low threshold parameters are tuned to the 
TDC hit efficiencies rather than trigger efficiencies 
as discussed above. 

• An alternate algorithm for simulating the PMT 
pulse shape was developed and tuned to observed 
PMT pulse characteristics. 

These differences manifest themselves as slightly differ- 
ent v e efficiency predictions and v c candidate rates in 
the data. Tests with radioactive sources have been per- 
formed to evaluate the quality of the second data recon- 
struction. The 22 Na and Am-Be efficiency runs shown in 
Fig. [ll] were reconstructed by the second analysis to test 
its efficiency prediction throughout the detector. The ra- 
tio of predicted to observed efficiencies over all the e + 
and neutron runs for the first reconstruction (1) and the 
second reconstruction (2) are plotted in Fig. [l4| While 
the results presented in the analyses below come from 
the first reconstruction code described above, the devel- 
opment of a second simulation and event reconstruction 
offers a useful crosscheck of the systematic uncertainties 
of the results. The differences between the two analy- 
ses were used to corroborate the estimate of systematic 
errors. 

IV. v c SELECTIONS AND BACKGROUNDS 

A. v a selection 

The trigger rate for time-correlated events (two triples 
occurring within 450 /xs) is ~1 Hz. Most of those events 
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FIG. 14. The ratio of predicted to observed efficiency of 
22 Na and Am-Be for the two reconstruction codes (1) and 
(2). 



are random coincidences of two uncorrelated triple hits, 
which occur individually at a rate of ^50 Hz, mostly from 
natural radioactivity. In order to select neutrino events, 
the following offline cuts are applied: 

• The energy reconstructed in both prompt and de- 
layed triples has at least one hit with E >1 MeV 
and at least two additional hits with E >30 keV. 
No single hit was allowed to be greater than 8 MeV. 

• The prompt triple is required to resemble a 
positron, i.e. annihilation j's each less than 600 
keV, and together less than 1.2 MeV. (This cut is 
the only one which treats the two triples asymmet- 
rically) . 

• At least one of the two triples in the event has more 
than 3.5 MeV of reconstructed energy for rejection 
of 7 backgrounds. 

• The prompt and delayed portions of the event are 
correlated in space and time (within 3 columns, 2 
rows, one meter longitudinally, and 200 /is). 

• The event started at least 150 fis (~5 neutron cap- 
ture times) after the previous veto tagged muon 
activity. 

The trigger and selection efficiencies are summarized in 
the first two columns of Table || 

In addition to corrections for selection cut efficiency 
and trigger efficiency, detector livetime is a substantial 
correction to the number of neutrinos seen and deserves 





Cumulative 




Data 


Cut 


efficiency 


Rate (d^ 1 ) 




1998 


1999 


1998 


1999 


Trigger 


0.271 


0.328 


69k 


106k 


Selection Cuts 


0.149 


0.177 


lk 


1.2k 


Live M vcto 


0.102 


0.121 


<50, see 


Live D aq- 


0.075 


0.112 


Table [v| 



TABLE I. Summary of efficiency and data rates after trig- 
ger, selections, and live time. The entry "Live MV eto" refers to 
the dead-time induced by the 150 /xs blanking window that 
accompanies each muon detected in the veto counter. 



some comment. Deadtime comes from two sources, the 
DAQ and the muon veto. DAQ livetime is the ratio of 
the number of triples the DAQ was available for digitizing 
to the total number of triples the trigger saw. These 
numbers are available from trigger scalers. The trigger 
livetime was measured to be >99.9%. The DAQ live time 
varies with the triple rate, and for the four data periods 
was determined to be 73.2%, 74.4%, 92.3%, and 91.8% 
for 1998 full power, 1998 refueling, 1999 full power, and 
1999 refueling, respectively. The higher livetime in 1999 
is due to improvements made in the trigger conditions. 

The muon deadtime can be further divided into two 
contributions: 150 /is of deadtime caused by each muon, 
which at = 1990 Hz left the detector live 74.2% of 
the time; and muons which interrupted a neutrino event 
between the positron and the neutron capture, which es- 
timated from the fit parameters of the Monte Carlo cap- 
ture time left 92.5% of events uninterrupted. The total 
uncertainty in the calculation of detector deadtime is less 
than 1% 



B. Backgrounds 

Backgrounds can be separated into two types: corre- 
lated and uncorrelated. Uncorrelated background events 
are due to unrelated triple hits which randomly coincided 
in the time window allowed. Although most of the events 
collected were random coincidences, almost all of this 
type of background is removed by requiring at least one 
subevent to have more than 3.5 MeV of reconstructed 
energy. These events do not have a time correlation be- 
tween prompt and delayed subevents (inter-event time) 
characteristic of neutron capture. They have instead 
a longer time correlation determined by the probability 
that the veto detected no muon between the prompt and 
delayed random triples. At a 2 kHz muon rate, this back- 
ground is seen as a 500 /is tail under the normal neutron 
capture distribution. By looking at the inter-event times 
of the candidate v c events at longer time scales, this back- 
ground can be measured. 

The inter- event time distribution after all neutrino se- 
lections (except the time correlation cut) is shown in 
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Fig. 15. The Monte Carlo for a pure neutron capture tical errors. 



sample is empirically fitted to the sum of two exponen- 
tials. There are two time constants due to the inhomo- 
geneity of the target: neutrons which remain in the scin- 
tillator have a 27 [is capture time, whereas those which 
enter the acrylic have a longer capture time due to the 
absence of Gd. The data inter-event time distribution is 
fitted to a function of three exponentials with fixed time 
constants consisting of the Monte Carlo fit r's multiplied 
by a third time constant of 500 /is. Integrating the result- 
ing 500 /is exponential of the uncorrelated background in 
the signal region gives an estimate of 4.1±0.2 events per 
day, or 9% of the v e candidates being uncorrelated back- 
ground events. 
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FIG. 15. The time elapsed between the prompt and de- 
layed portion of neutrino candidate events for Monte Carlo 
simulation and data. The Monte Carlo is fit to two exponen- 
tials. The data are fit to three exponentials, constrained to 
have the time structure of the simulation and an additional 
contribution due to uncorrelated background (dashed line). 

To measure the uncorrelated background in smaller 
parts of the data set, the statistical accuracy of the three 
exponential fit method becomes unacceptably poor. A 
simpler method is therefore used in conjunction with the 
above fit. For inter- event times longer than 200 [is, the D e 
candidates are dominated by uncorrelated backgrounds. 
The integrated number of candidates from 200-400 /is 
is scaled to estimate the number underneath the signal 
region (<200 /is). Using the scaling from the fit of the 
entire data set shown in Fig. [l5|, the uncorrelated back- 
ground was measured in approximately month-long in- 
tervals as shown in Fig. [l6| For both the 1998 and 1999 
data sets, the rates are found to be stable within statis- 
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FIG. 16. The uncorrelated background rate as measured in 
approximately month-long segments of the data. Overlain is 
a fit to a constant within each year's data; the fit quality is 
consistent with a stable background. 

Correlated backgrounds have the neutron capture 
inter-event time structure of the v e candidates. These 
events come mainly from cosmic muon induced fast 
neutrons from spallation or muon capture, as shown 
schematically in Fig. [Tt]. These fast neutrons can either 
(1) induce more neutrons via spallation, two of which 
can be captured in the detector with one capture mim- 
icking a positron signature; or (2) they can cause pro- 
ton recoil patterns in the central detector which appear 
as a positron signature and then get captured. Spalla- 
tion neutrons originate from muons passing through the 
walls of the lab without hitting the veto detector or from 
muons passing through the detector shielding undetected 
by the veto. Muon capture neutrons mainly originate 
from muons stopping in the water buffer without regis- 
tering in the veto. 

To illustrate some properties of correlated background, 
Fig. |l8| shows the time elapsed since the previous veto hit 
for v c candidates, with all selection cuts applied except 
that of the previous muon timing. This distribution is fit 
to a three exponential function analogous to that used for 
the inter- event time fits. The two time constants for neu- 
tron capture are not identical to those for v c events, but 
tend to be smaller since after passage of a muon there 
are often more than one neutron in the detector to be 
captured. The third exponential time constant is again 
constrained to 500 /is as expected in a random sampling 
of events unrelated to the previous muon. Since at very 
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FIG. 17. A schematic drawing of the detector showing a signal event and two examples of interactions which contribute to 
the correlated background. Neutrons which cause correlated backgrounds come mainly from muons undetected by the veto 
and captured in the water buffer (left), and spallation from muons in the lab walls (right). These neutrons in turn can either 
induce more neutrons via spallation (as shown at right) or cause three coincident proton recoils as they thermalize (as shown 
at left). 
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FIG. 18. The time between the previous muon traversing 
the veto counter and a neutrino candidate event trigger. The 
fit is to a sum of three exponentials. The dashed portion 
denotes events uncorrelated to the previous muon, the rest are 
correlated background events induced by cosmic muons. The 
neutron capture time constants, shorter than the simulation 
prediction for v e events, imply multiple neutron production. 

short times there are other contributions such as muon 
decay, times less than 15 /is are excluded from the fit. 



Muon-induced-neutron backgrounds dominate the candi- 
dates in the first 150 /xs after the previous tagged muons, 
motivating the selection cut on \i timing. 

In order to show that the correlated background was 
constant in time, the previous muon time cut was dis- 
abled and a plot was made of the v candidate rate ver- 
sus time, as shown in Fig. [l9|. When fit to a constant for 
each year, a x 2 /n.d.f. of 382.7/371 is obtained, which has 
a 33% likelihood, indicating that the detector efficiency 
for correlated backgrounds was stable during each year's 
data taking. 

Aside from the detector efficiency for background, how- 
ever, a loss of veto efficiency could also cause a fluctua- 
tion in background. (The rates fit in Fig. |l^ are with the 
muon timing selection disabled, and hence do not vary 
with veto inefficiency.) To track veto efficiency, the veto 
hit patterns recorded with each event are used. If a muon 
hit was recorded only on the bottom of the veto, where 
only exiting muons are seen, then the muon must have 
entered the veto without recording a hit. By measuring 
the percentage of these events a one hit missed veto in- 
efficiency of (4 ± 1)% is found as mentioned above. The 
through-going [two hits missed) veto inefficiency is mea- 
sured to be 0.07±0.02% by looking at the rate of /i tracks 
triggered on in the central detector. These inefficiencies 
were tracked in time to assure their stability. 

C. Neutron— i> c direction correlation 

The neutrons produced in the inverse beta decays 
will have momenta slightly biased away from the source, 
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FIG. 19. The number of candidates per day in the 373 
runs used in this analysis, fit to a constant for each year. 
All neutrino selections except the previous muon time cut are 
used, leaving a data set dominated by muon-induced-neutron 
background. The different rates in 1998 and 1999 are due to 
trigger and livetime improvements. 

whereas no correlation is expected for background. This 
effect is the consequence of momentum conservation 
which requires that the neutron should always be emitted 
in the forward hemisphere with respect to the incoming 
D c . Such a correlation has been observed already in the 
Gosgen experiment |||] and again at Chooz The 
theoretical treatment of the effect can be found in |L7| . 

The signal to background ratio can be independently 
verified using this effect. The v c source is to the left of 
the detector in Fig. |l|. The relative horizontal location 
(relative column in the target cell array) of neutron cap- 
ture cascade cores versus positron ionizations for data 
and the simulation of the v c signal are plotted in Fig. |2C| . 
Defining the asymmetry Adata = fef; in terms of the 
number of neutrons captured one column away from the 
source R and one column toward the source L, a slight 
asymmetry 0.050 ± 0.017 is found in the data, at 2.9 a 
significance. Using the Monte Carlo simulation which 
gives Amg = 0.134 to estimate the portion of the data 
consisting of P e signal and assuming the background to 
be symmetric in this variable, an effective signal to noise 
ratio 
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is found. This value agrees well with the ratio of 
0.81 ± 0.03 found with the swap analysis method de- 
scribed below. 
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FIG. 20. The relative horizontal location of the neutron 
capture and positron ionization in the detector target. The 
kinematics of the inverse beta decay cause a slight bias in the 
momentum of the neutron from inverse beta decay, seen here 
as an asymmetry in the populations of neutrons captured one 
column away or towards the v e source. Backgrounds should 
not exhibit such asymmetry. The Monte Carlo prediction for 
the pure u c signal, normalized to unity, is also plotted. 



V. ANALYSIS 

The data set presented here was taken from July 1998 
to September 1999 in 373 short runs, each on average 
about 12 hours long. In 1998, 35.97 days of data were 
taken with the three reactors at full power and 31.35 days 
with one of the reactors at a distance of 890 m off for re- 
fueling. The detector was then taken offline in Jan/Feb 
of 1999, when DAQ improvements were made to increase 
livetime, and the high trigger thresholds were lowered by 
30% to increase trigger efficiency The 1999 data set in- 
cludes 110.95 days with all three reactors at full power 
and 23.40 days with the 750 m baseline reactor off for 
refueling. Thus, the entire data set has four distinctly 
different periods, with three different baseline combina- 
tions and neutrino fluxes. 

After all selection cuts there is still substantial back- 
ground in the remaining data set. The correlated back- 
ground, coming mainly from muon induced neutrons, is 
difficult to predict and subtract. The yield and spectrum 
of neutron spallation is a function of muon flux and en- 
ergy, which in turn is a function of depth. While some 
measurements of fast neutron spectra and fluxes have 
been done in the past, there is no model which can con- 
sistently predict the fast neutron production. Below we 
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present three methods used to extract the v c signal from 
data. 



the uncorrelated background was measured and removed 
from the data before the subtraction. 



A. Analysis with the on-off method 

The conceptually simplest method of subtracting back- 
ground is to take advantage of periods of reduced power 
levels of the reactor source. Ideally all three reactors 
would be down at once allowing for a direct measurement 
of the background. However, in practice only one of the 
three Palo Verde reactors was refueled at any given time. 
These reduced power periods occurred twice annually for 
about a month. Each year's data set is treated indepen- 
dently, subtracting 1998 off from 1998 on and 1999 off 
from 1999 on, since the efficiency of the detector changed 
between the two years. By subtracting these data taken 
at reduced flux from the full flux data, a pure neutrino 
sample is retrieved albeit containing the statistical power 
of only a small portion of the potential data set: the sub- 
traction is limited by the two months of refueling time 
and treats the v e flux from the two reactors still at full 
power as background. 

The primary concern arising from use of this method, 
aside from the loss of statistics, is guaranteeing that the 
background rates during the on and off periods were sta- 
ble. Both correlated and uncorrelated backgrounds were 
carefully tracked to ensure stability as discussed above. 





1998 


1999 


Mm) 


890 


750 


ON iV cand (day- 1 ) 


38.2 ± 1.0 


52.9 ±0.7 


OFF iV cand (day- 1 ) 


32.2 ± 1.0 


43.9 ± 1.4 


ON-OFF iV cand (day" 1 ) 


6.0 ± 1.4 


9.0 ±1.6 


Total efficiency ON (OFF) 


0.0746 (0.0772) 


0.112 (0.111) 


Robs (day -1 ) 


95 ±19 


77 ± 14 


7? ca i c (day -1 ) 


63 


88 



TABLE II. Results for the simple on-off background sub- 
traction analysis, showing candidate rates in the data N and 
efficiency corrected v c interaction rates R observed and cal- 
culated. The data sets for each year are considered indepen- 
dently here. Uncertainties are statistical only. Systematic 
errors are estimated to be 10%. 

The numerical results of this analysis of the total rate 
are summarized in Table ||. After correcting for efficiency 
(for the no-oscillations scenario) and livetime, the data 
sets were subtracted to find observed neutrino interaction 
rates in the detector. No significant deviation from the 
expected neutrino interaction rates was found at either 
baseline distance. 

The results from the alternate reconstruction (2) for 



this analysis are shown in Table III for comparison. This 
analysis selects about 5% more candidates, but also gives 
a correspondingly higher efficiency. For this analysis, 





1998 


1999 


L(m) 


890 


750 


ON 7V cand (day" 1 ) 


37.3 ± 1.2 


49.3 ±0.7 


OFF 7V cand (day" 1 ) 


31.6 ± 1.2 


38.6 ± 1.6 


ON-OFF iV cand (day" 1 ) 


5.7 ± 1.7 


10.7 ± 1.8 


Total efficiency ON (OFF) 


0.0809 (0.0838) 


0.121 (0.121) 


R obs (day -1 ) 


85 ±20 


89 ±15 


#caic (day -1 ) 


63 


88 



TABLE III. Results from the alternate reconstruction (2) 
for the simple on-off background subtraction analysis. Uncor- 
related background was accounted for in the the data before 
subtraction. Uncertainties are statistical only. Systematic 
errors are estimated to be 10%. 



In order to test the results for oscillation hypotheses in 

,2 „;„2, 



the two flavor Am 2 — sin 2 28 plane, a x 2 analysis is per- 
formed comparing the calculated R C a\c,ij an d observed 
Robs.ij spectra divided into 1 MeV bins j for each year 
i. The spectra used are the prompt energies of the two 
subtracted data sets. At each point in the oscillation pa- 
rameter plane, taking into account the changes in detec- 
tor efficiency due to distortions of the neutrino spectrum, 
the quantity 



^2 _ ( a Rcnlc,ij — R~obs,ij) 2 _j_ (a — l) 2 



i=l 



a: 



(5) 



i j 



syst 



is computed, where a accounts for possible global nor- 
malization effects due to systematic uncertainties (dis- 
cussed below) across both periods and ay is the statis- 
tical uncertainty in each bin. Systematics which can af- 
fect spectral shape, mainly energy scale uncertainty, are 
negligible relative to the statistical uncertainties in the 
analysis. The function is minimized with respect to a. 
The point in the physically allowed parameter space with 
the smallest chi-square Xbcst was found, which represents 
the oscillation scenario best fit by the data. 

The 90% confidence level (CL) acceptance region is 
defined according to the procedure suggested by Feldman 
and Cousins |2^] by: 

AX 2 - X 2 (Am 2 , sin 2 28) - X Lst > Ax^A™ 2 , sin 2 20) 

(6) 

where x 2 (Am 2 , sin 2 28) is the minimized fit quality at 
the current point in Am 2 — sin 2 28 space and Ax 2 r j t is 
the CL x 2 cutoff. Due to the sinusoidal dependence of 
the expected rates on the oscillation parameters and the 
presence of physically allowed boundaries to those pa- 
rameters, the cutoff is not simply the Ax 2 one would 
analytically find for a three parameter minimization but 
has to be calculated for each point in the plane. To find 
the Ax 2 r i t for a point, the experiment is simulated 10 
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times under the assumption that the oscillation hypothe- 
sis represented by that point is true. For each simulated 
data set, a xtest is extracted and a Ax 2 found for the 
point. These 10 4 Ax 2 , the simulations' fit qualities to 
the hypothesis, are then ordered. The Ax 2 of which 90% 



to a very high degree of confidence by the independent 
efficiency calibrations of the detector discussed in pre- 
vious sections. In addition this result has no effect on 
the exclusion plot in Fig. ^l] because, as shown in Ta- 
ble VI the no-oscillation hypothesis has actually better 
of the simulations are a better fit is a 90% CL and there- x 2 / n -d-f- than the minimum. Also the exclusion plots 



fore that oscillation hypothesis' Ax 2 r j t . 

The region excluded by the analysis is shown in curve 
(a) of Fig. EH]. The results of the x 2 analysis, including 
the oscillation parameters' best fit to the data, are sum- 
marized in the first column of Table VI further below. For 
the on- off analysis a best fit preferring the no-oscillation 
hypothesis was found. 

In addition to the analysis of the absolute P e rates ob- 
served, one can analyze the shape of the spectrum of neu- 
trinos seen independently of the absolute normalization, 
thereby relieving the result of most systematic uncertain- 
ties. The x 2 is calculated at each point in the oscillation 
parameter plane as in Eqn. ( |5|) , with no constraint on nor- 
malization (<7 S yst — > oo). The same procedure as before is 
followed in defining a 90% CL region in the sin 2 20-Am 2 
plane. At large Am 2 where T> e of all energies are oscillat- 
ing many times within the baseline, the energy spectrum 
of the incident flux is affected only in magnitude. As a 
result, the region excluded in the plane does not extend 
to large Am 2 , as shown in Fig. Efl (b). 



> 
u 



10 



10 



10 



10 



.1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 
: (a); 


1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 . 

(b) --: 






■ ... 1 .... 1 .... 1 .... 1 .... 1 





0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

sin 2 29 

FIG. 21. 90% CL limits of the on-off analysis, (a) on-off 
subtraction with normalization constraint, and (b) using only 
the spectral shape. 

For the spectrum analysis, when the normalization a 
is left free, the minimum x 2 is obtained for a ~ 2 (see 
Table VI) and maximum mixing. This is clearly an un- 



bascd on Eqn. (|6|) and either Xbcst or Xno.osc. are found 
to be virtually identical. Furthermore, changing the bin 
size from 1 MeV to 0.5 MeV does not appreciably change 
the exclusion plot, either. 

Since the analyses reported above and in the follow- 
ing sections finds no evidence for neutrino oscillations, 
the spectra of the two years are added and the summed 
spectrum is plotted in Fig. along with the Monte Carlo 
expectation. 



100 
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5 6 7 
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FIG. 22. Total energy of the prompt subevent (~ E„ a — 1.8 
MeV) in 1 MeV wide energy bins (observed and expected for 
no oscillation) from the on-off analysis, summed over 1998 
and 1999. The spectra are in agreement with each other, 
with a x 2 / n -d.f. = 8.5/8. Uncertainties are statistical only. 



B. Analysis with the reactor power method 

Part of the statistical limitations of the direct subtrac- 
tion of the preceding analysis is a result of the separa- 
tion of the data set by year. By correcting the four pe- 
riods for efficiency and then subtracting the respective 
reduced flux from full flux periods, the subtraction is 
forced to treat the P e flux from two of the reactors as 
background. A second x 2 analysis was performed which 
effectively uses the full v e flux of the refueling periods. 

To use the 1998 and 1999 data sets together, the change 
of both signal and background efficiencies are accounted 



physical result since such large value of a can be excluded 
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Events 


1998 ON 


1998 OFF 


1999 ON 


1999 OFF 


(day" 1 ) 




(890 m) 




(750 m) 




38.2±1.0 


32.2±1.0 


52.9±0.7 


43.9±1.4 


b 


19.5 ± 1.7 


26.3 ±2.2 1 


Adctected 


18.7±2.0 


12.7±2.0 


26.6±2.3 


17.6±2.6 


Robs 


225±24 


140±22 


216±19 


140±21 


Rcalc 


218 


155 


218 


130 



TABLE IV. The rate of candidates iV can d, signal iVdetected, 
and background b found by the reactor power analysis, us- 
ing a \ 2 analysis of the data set run by run. The effi- 
ciency-corrected total interaction rates R are also listed. Er- 
rors are statistical only. The systematic uncertainty is esti- 
mated to be 10%. 



for. The v c efficiency difference is found through the de- 
tector Monte Carlo simulation. The efficiency change for 
background from 1998 to 1999, which is not necessarily 
the same as for P c , is extracted via the high statistics 
fits of correlated background shown in Fig. [TJ. The 27% 
increase in background efficiency observed roughly corre- 
sponds to what the Monte Carlo simulation predicts for 
a background composed mainly of double neutron cap- 
tures. Uncorrelated background accounts for less than 
10% of the candidate set; this background's efficiency 
changed by a similar amount within the measurement 
statistics seen in Fig. [jl]. 

The combined data sets are analyzed for oscillation 
hypotheses by calculating the x 2 summed over runs i 
(runs with less than ten candidates are combined with 
adjacent runs): 



+ b)-N cand , l f , {a -If 



(7) 



y syst 



where iV can d,i is the total v c candidate rate, N c& \ c ,i is 
the calculated rate, a is the overall normalization as be- 
fore, and b is the background rate. The background, b, 
is scaled as appropriate for the year but is otherwise as- 
sumed to be constant. The function is minimized at each 
point with respect to b and a. We found no evidence 
for oscillations and the 90% CL plot, shown in Fig. |2^, 
curve (a), is then constructed around the Xbest as before 
by comparing A\ 2 with Axl rit at each point. The pre- 
dictions of signal and background from this fit for the 



no-oscillation hypothesis are shown in Table IV. The no- 
oscillation likelihood and best fit results of the x 2 analysis 
with this method are summarized in the third column of 
Table [v| 

C. Analysis with the swap method 

A third analysis is used which has the potential of using 
the full statistical power of the neutrino data set by sub- 
tracting background directly. The method, discussed in 



more detail elsewhere p6| , takes advantage of the asym- 
metry of the prompt (positron) and delayed (neutron 
capture) subevents of the neutrino signal. The data se- 
lection and trigger treat the two portions of the event 
identically with the exception of two cuts designed to 
isolate events with annihilation- like 7's in the prompt 
triple. 

The candidates remaining after the selection cuts can 
be written as: 



N = B„ 



D 



B 



pn 



S u (8) 

where B unc , B nn , and B pn are uncorrelated, two- neutron, 
and proton-recoil-neutron-capture backgrounds respec- 
tively; and S v is the neutrino signal. Applying the same 
neutrino cuts with the positron cuts reversed, or swapped, 
(such that the positron cuts are now applied to the de- 
layed triple) gives: 



N' = B„ 



B n 



e\B pn + ezSv 



(9) 



Since the uncorrelated background and two neutron cap- 
ture backgrounds are symmetric under exchange of the 
prompt and delayed triples, their efficiencies with the re- 
versed cuts remain the same. The parameters t\ and £2 
denote the relative efficiency change for proton recoils 
and neutrino signal under the swap, respectively. 

The positron cuts are highly efficient for positron anni- 
hilation events but have poor efficiency for neutron cap- 
tures. The Monte Carlo simulation is used to estimate 
£2 = 0.159. Subtracting (@) from ([)[) leaves the majority 
of the neutrino candidates and only proton recoil back- 
ground: 



N - N' = (1 - e 1 )B pn + (1 - e 2 )S v . 



(10) 



To estimate (1 — e\)B pn , it is noted that the proton 
recoil spectrum extends beyond 10 MeV, well above the 
positron energies of the neutrino signal and other sources 
of background. These measured high energy events can 
be used to normalize the B pn background in the signal 
using the Monte Carlo ratio: 



= flMC (Ei|C+ <8Me V) 
r -fl«C(E li6+ > lOMeV)' 



(11) 



where 



(E l e + < 8MeV) is the fraction of simu- 
lated B pn events passing the normal v c selections, and 
B p v ^ c (E l o + > lOMeV) are the fraction of simulated 
events in the high energy background region. Multiply- 
ing the ratio r by the measured high energy proton recoil 
rate gives the B pn background contribution: 



B pD = rB^ ta (E liC+ > lOMeV). 



(12) 



The neutrons which cause the proton recoil back- 
ground are created either by muon capture or spallation 
in the laboratory walls, or by muons entering the veto 



1G 



Period 


1998 ON 


1998 OFF 
890 m reactor off 


1999 ON 


1999 OFF 
750 m reactor off 


time (days) 


35.97 


31.35 


110.95 


23.40 


i> e overall efficiency (%) 


7.46 


7.72 


11.2 


11.1 


£C a (E liC+ > 10MeV)(day- 1 ) 
(1 — ei)B pll (day -1 ) fi spallation 


8.79 


9.09 


13.52 


13.29 


-0.88 


-0.91 


-1.35 


-1.33 


(1 — ei)B pn (day -1 ) n capture 


0.58 


0.58 


0.86 


0.86 


N (day" 1 ) 


38.2 ± 1.0 


32.2 ±1.0 


52.9 ±0.7 


43.9 ± 1.4 


N' (day -1 ) 


24.6 ± 0.8 


21.2 ± 0.8 


32.3 ± 0.5 


31.7 ± 1.2 


N„ (day" 1 ) 


16.5 ± 1.4 


13.5 ±1.4 


25.1 ±0.9 


15.0 ±1.9 


Total background B unc + B nn + B pn (day x ) 


21.7 ± 1.0 


18.7 ±1.0 


27.8 ±0.6 


28.8 ± 1.3 


Robs (day -1 ) 


221 ± 19 


174 ± 17 


225 ± 8 


137 ± 17 


-Rcaic (day -1 ) 


218 


155 


218 


130 



TABLE V. Results for the swap analysis, including the various background estimates. Uncertainties are statistical only. 



counter undetected. The spectrum of the fast neutrons 
from spallation is not well understood. However, such 
spectrum can be decoupled somewhat from the resulting 
proton recoil spectrum. The expected backgrounds were 
simulated for various possible fast neutron spectra and 
the resulting e\ and r for neutrons created in the lab walls 
were calculated. The same calculation was performed for 
neutrons created in the passive detector shielding by un- 
tagged muons; in this case, the expected yield is much 
smaller, being only a few percent of that from the walls. 
The simulated spectra of spallation neutrons are chosen 
to span the wide range of predictions quoted in literature. 

A value for e\ of 1.14 ± 0.07 is found after averaging 
over spectra, implying that the spallation proton recoil 
background is essentially symmetric like the other back- 
grounds. Upon simulating the possible spectra, the quan- 
tity (1 — e\)r — 0.1 ± 0.05 is found to vary little. 

The yield and spectrum of neutrons from muon cap- 
ture are reasonably well understood. Since these neu- 
trons tend to be lower in energy, only those created in 
the vicinity of the detector have any efficiency for creating 
background. Knowing the veto inefficiency to miss stop- 
ping muons (4±1)%, the capture rate in water surround- 
ing the detector and its contribution to the background 
can be estimated using Monte Carlo simulation. Overall 
this proton recoil background appears to be symmetric 
as well, ei = 0.77 ± 0.32, meaning that the subtraction 
also strongly rejects this background. The uncertainty 
of the residual background (1 — e\)B pn is conservatively 
estimated to be about 160%, corresponding to ~4% error 
on iV Po . 

The results of this analysis are summarized in Table 
Overall 



analysis but uses the background estimated by the swap 
method subtraction instead of minimizing the function 
with respect to background. Again, we find no evidence 
for neutrino oscillations and the excluded region for this 
analysis is shown in Fig. p3|, curve (b). 
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10 




Rcalc 



1.04 ± 0.03(stat.) ± 0.08(sys. 



(13) 
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sin 2 26 

FIG. 23. 90% CL limits of oscillation parameters for (a) 
the reactor power analysis, which fits for background based on 
varying power levels, and (b) the swap analysis, which directly 
subtracts background from the data set. The Kamiokande 
v B -Vp, atmospheric neutrino 90% CL (dashed line) and best 
fit (star) are also shown. 



D. Systematic uncertainties 



The background estimates returned by the reactor power 
analysis in Table |iy| compare well with the results of the 
swap analysis. The 90% CL region for this analysis fol- 
lows the same \ 2 formula, Eqn. m, as for the reactor power 



The systematic uncertainties have three sources: the 
prediction of expected P e interactions, the efficiency esti- 
mate, and, for the swap analysis, the B pn estimate. The 
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Analysis 


On-Off 


Spectrum 


Power 


Swap 


Xbcst/n.d.f. 


17.9/13 


16.9/13 


317.6/325 


0.9/1 


Qbcst 


1.00 


1.99 


1.02 


1.03 


Xno.osc/n.d.f. 


17.9/15 


17.9/15 


317.7/327 


0.9/3 



TABLE VI. The best fit to the data for the oscillation pa- 
rameters of the two flavor model for each of the analyses pre- 
sented here: on-off, spectral shape, reactor power, and swap. 
The x 2 °f the no-oscillation hypothesis is also tabulated for 
each analysis. For the reactor power analysis all runs with 
less than 10 u c candidates were combined with adjacent runs, 
leaving 327 d.o.f. See text above for an explanation of the 
"Spectrum" analysis. 



expected v e uncertainty is dominated by the conversion 
of fission rates into neutrino fluxes, which relies on direct 
empirical measurements of (3 spectra emitted by the iso- 
topes. The Bugey experiment p7[ , which directly mea- 
sured the neutrino flux and energy spectrum emitted by 
a reactor at short baseline, found agreement within 3% 
using the same methods; the 3% value is used here as the 
estimated uncertainty. 

The efficiency uncertainty can be further subdivided 
into that arising from direct comparisons of Monte Carlo 
e + and neutron efficiency from calibration measurements 
and that arising from the selection cuts themselves. The 
calibration runs taken with the positron and neutron 
sources, when compared with Monte Carlo simulations, 
shows overall agreement across all locations of better 
than 1% in the efficiency predictions. However, the run- 
by-run agreement was at a level of 4% for neutrons and 
11% for positrons. Since the 22 Na source is similar to 
the inverse beta decay signal with the e + close to de- 
tector threshold, the positron efficiency uncertainty over 
the entire v c spectrum was estimated to be closer to 4% 
in any particular location. These run-by-run variations 
are then used as our systematic uncertainties in the v e 
efficiency. 

To test the robustness of the event selection, each cut 
is varied within a reasonable range and variations of the 
ratio between data and Monte Carlo are examined. In 
order to take into account correlations all cuts were var- 
ied simultaneously by randomly sampling a multidimen- 
sional cut space. The rms of the resulting ratio of ob- 
served/expected is given as the selection cut uncertainty. 

The swap method analysis has a somewhat smaller un- 
certainty for the selection cuts variation as the subtrac- 
tion tends to cancel out systematics. However, the swap 
analysis uses a Monte Carlo estimate of the proton recoil 
background. Due to limited Monte-Carlo statistics and 
the uncertainty in the fast-neutron energy spectrum, a 
4% uncertainty is assigned to the neutrino signal. All 
of the systematic uncertainties are summarized in Ta- 
ble VII. The total systematic uncertainty is obtained by 
adding the individual errors in quadrature. 

The development of a second simulation and event re- 



Error Source 


On Minus Off(%) 


Swap(%) 


e + efficiency 


4 


4 


n pffipipTiPv 


3 


3 


P e flux prediction 


3 


3 


v e selection cuts 


8 


4 


S pn estimate 




4 


Total 


10 


8 



TABLE VII. Summary of the systematic uncertainties. 



construction proved to be helpful in understanding sys- 
tematic uncertainties of the analyses due to the algo- 
rithms chosen. For comparison the results for the on-off 
analysis from both reconstructions are shown in Tables |FJ 



and III . An independent analysis of systematic errors was 
performed for the second reconstruction, similar to the 
method described above, giving comparable results. 



VI. CONCLUSION 

In conclusion, the data taken thus far from the Palo 
Verde experiment show no evidence for D c — > v x oscil- 
lations. This result, along with the results reported by 
Chooz and Super-Kamiokande excludes two fam- 
ily v^-Vc mixing as being responsible for the atmospheric 
neutrino anomaly as originally reported by Kamiokande 
jD. Later results of Super-Kamiokande, in particular 
data on the zenith angle distribution of muons and elec- 
trons, suggest that muon neutrinos strongly mix with 
either v T or with a fourth flavor of neutrino sterile to weak 
interaction. Clearly it is becoming important to include 
at least three neutrino flavors when studying results from 
oscillations experiments. 

The most general approach would involve five un- 
known parameters, three mixing angles and two inde- 
pendent mass differences. However, an intermediate ap- 
proach consists of a simple generalization of the two fla- 
vor scenario, assuming that m\ 3> m 2 ,™ 2 . (i.e. Ar 

2 _ TTT v.;l^ A™2 



'13 



A 



'23 



while Ar 



(, 2 2 ~ 0). In such a case the mix- 
ing angle Q\i becomes irrelevant and one is left with only 
three unknown quantities: Am 2 , #13 , and 6*23 . With this 
parameterization the v e disappearance is governed by 



P{v e — > v x ) — sin 2 26>i3 sin' 



Am 2 L 
4£„ 



(14) 



while the Vp, — » v T oscillations in this scenario responsible 
for the atmospheric neutrino results, are described by 



Am 2 I j 

v T ) = cos 4 sin 2 2#23 sin 2 . (15) 



A preliminary analysis of the atmospheric neutrino 
data based on these assumptions has been performed pgj ] 
and its results are shown in Fig. ^ for the v e disap- 
pearance channel. One can see that while the relevant 
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Hungarian OTKA fund, and another (L.M.) from the 
ARCS Foundation. 
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sin 2 29 13 

FIG. 24. Exclusion plot showing the allowed region of 9is 
and Am 2 based on the Super-Kamiokande preliminary analy- 
sis (the region inside the dotted curve). The region excluded 
by the neutrino reactor experiments are to the right of the 
corresponding dashed and continuous curves. 



region of the mass difference Am 2 is determined by the 
atmospheric neutrino data, the mixing angle $13 is not 
constrained very much. Here the reactor neutrino results 
play a decisive role. 

We plan to continue taking data through Summer of 
2000, which will provide two additional reduced flux re- 
fueling periods. 
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